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The analytical results obtained in the infinite mass and strong coupling limits of QCD are difficult to reconcile 
with the predictions of the Monomer Dimer Polymer algorithm. We have reconsidered in detail the results 
obtained with this simulation scheme and evidences of severe convergence problems are presented for the SU{3) 
and 5(7(2) gauge group cases. 
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Introduction 

Finite Density QCD is affected by the well 
known sign problem that has prevented, in most 
cases, any success in simulate this theory and, 
until now, no solution is at sight. The only ex- 
ceptions in this scenario are strong coupling simu- 
lations performed with the Monomer Dimer Poly- 
mer (MDP) algorithm ||l|. This algorithm is able 
to provide results not affected by early onset and 
in partial agreement with the Mean Field (MF) 
prediction. More recently another theoretical ad- 
vance has been achieved solving exactly the P — 
QCD in the limit of infinite mass and chemical 
potential In the next section we discuss the 
compatibility of MDP results with the infinite 
mass SU{3) solution; indeed we have found dif- 
ficult to conciliate the numerical and analytical 
predictions. This leds us to reconsider the MDP 
algorithm more carefully and we have found evi- 
dences of convergence problems. 

We have also considered the SU{2) case where 
the sign problem is not present. Results obtained 
using the Gran Canonical Partition Function 
(GCPF) formalism turn out to be in very good 
agreement with the Hybrid Montecarlo (HMC) 
calculations |^ while, once again, MDP results 
Q are inaccurate in the critical region. 



SU{3) case 

The main feature of infinite mass QCD I IS 
that at zero temperature (it — + oo) the system 
undergoes a first order saturation transition and 
the phase of the Dirac determinant is not rele- 
vant. However at non-zero temperature (Lt fi- 
nite) the system has only a smooth crossover and 
the phase is relevant. The MDP results are some- 
how puzzling if considered in the light of the in- 
finite mass solution. In [|l| Karsch and Miitter 
saw a strong first order transition for = 4 
and is = 4,8 with m G [0.1,0.7], while we know 
that for m — > oo the number density is a smooth 
function of fi for any finite Lj. To reconcile the 
m oo solution and MDP results we have to 
suppose that the transition disappears at some 
large bare mass m, or that the infinite mass limit 
of QCD is singular. These statements seem both 
unnatural. The former would imply the existence 
of a (large) physical scale where the system be- 
haviour changes, washing out the transition. The 
latter is unplausible too, because the Dirac deter- 
minant approaches its infinite mass limit contin- 
uously. In order to have data more easily compa- 
rable with the analytic predictions we have tried 
to use the MDP algorithm |^ directly in the large 
mass regime. The authors of the MDP code no- 
ticed M that for small masses (i.e. m < 0.1) the 
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Figure 1. Number density as function of /i with: 
fi = start (diamonds), saturated start (squares). 



algorithm becomes non-ergodic. When we tried 
MDP simulations for m > 1.0 we saw a similar 
behaviour: for no value of ^ the system moved 
into the saturated phase. Moreover a degrada- 
tion in performance (for all mass values) has been 
observed increasing the lattice size but, what is 
more surprising, this degradation seems to be re- 
lated only to the value of Lf. We have not been 
able to perform simulations for > 4. 

The behaviour of MDP code prevents any di- 
rect comparison between the large mass results 
and MDP ones. The limited applicability of 
MDP algorithm raises doubts on its general va- 
lidity. We have repeated the MDP simulations at 
TO = 0.1 in a 4"^ X 4 lattice (values used in the 
original paper 0|), using as initial configuration 
either n(^) = or n(/i) = 1 and 0(10*^) Mon- 
tecarlo steps for each value of ^. From fig. 1 it 
is evident a signal of hysteresis in the data; the 
fi — start undergoes a strong saturation tran- 
sition at ^ = 0.69 (pubhshed result ||l|), while 
the run with saturated start jumps in the zero 
density phase at /i = 0.58. This result should 
lead us to cautioun on the determination of the 
critical point and may well reconcile the MDP re- 
sults with MF predictions (the MF critical point 
= 0.61 lies inside the hysteresis). To be con- 
fident with MDP results we should observe sev- 
eral flip-flops in the Montecarlo history to con- 



clude that a clear two state signal is present. We 
have considered runs of up to O(IO^) configura- 
tions and we did not succeed to see any flip-flop 
for any fi inside the hysteresis region. Starting 
from a zero density configuration nothing hap- 
pens until we get close to /z = 0.69. At this point 
the system has some probability to go in the satu- 
rated phase. Once the system is in the saturated 
phase it never goes back. The same behaviour 
has been observed (near fi = 0.58) for the satu- 
rated start. Varying the quark mass only changes 
the hysteresis position unless we reach too small 
(to < 0.1) or too large (to > 1.0) values. From 
these numerical evidences we may conclude that 
the hysteresis behaviour of the system is inde- 
pendent on the statistics for any value of /u. The 
MDP code seems to have convergence problems 
in the most interesting region of fj,, independently 
on the temporal lattice extent and the quark mass 
value. 

SU{2) case 

Let us now address the MDP convergence in 
the case of SU{2) gauge group. The motiva- 
tion is twofold: firstly we wonder if the problems 
present in SU (3) are universal independently on 
the gauge group. Moreover SU{2) ofl^ers us the 
possibility of using conventional simulation algo- 
rithms. In this case, in fact, quarks and anti- 
quarks belong to the same (real) representation: 
the Dirac determinant is real and positive also 
for non-zero fi. The main published MDP re- 
sults in SU{2) concern the number density and 
the chiral condensate in a 4^^ x 4 and 8"^ x 4 lat- 
tice at TO = 0.2 g. Other results for SU{2) at 
f] — and /i / are obtained using the HMC 
algorithm in a 4^ x 4 at m = 0.2 In order 
to have results with /i varying continuously, we 
have used the GCPF scheme, performing simula- 
tions with Lt — 4 and Lg = 4, 6, 8, at the quark 
mass values to = 0.1, 0.2, 0.4. In fig. 2 we report 
the number density and the chiral condensate ob- 
tained in our simulations compared with HMC 
results appeared in fig. 4 of [P; error bars are 
reported at some values of /i. It is evident that 
our simulation reproduces the HMC results quite 
accurately. This agreement between GCPF and 
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Figure 2. Number density and chiral condensate 
as a function of /i in a 4^^ x 4 at to = 0.2 from 
GCPF (line) and HMC (symbols). 



Figure 3. Chiral condensate at to = 0.2 in 4^^ x 4 
(dots), 6^ X 4 (dashes), 8^ x 4 (line) and MDP 
results in a 4^ X 4 (diamonds) and 8^ x 4 (squares) . 



HMC suggests that sampling problems, found in 
SU{3) with GCPF (I, are not present in S'[/(2). 
Simplified models predict for SU{2), at least at 
small temperature, a phase transition at half of 
the mass of the lightest baryon of the theory (de- 
generated with the pion at n — 0). We have 
then computed the pion mass in a 6'^ x 12 at 
TO = 0.1,0.2,0.4. To extract the critical /i we 
have used the following critcrium. The number 
density appears, increasing the volume, almost 
zero up to /ic, with a linear rise beyond it and 
flat at large fi (saturation). To identify the crit- 
ical point we have computed dn{ii)/dfi for two 
different volumes and defined fic as the position 
of the first crossing of the curves. In the limit 
y — > oo this defines correctly the fi at which the 
linear behaviour begins. In table we report our 
fic and one half the pion mass for different values 
of to: 



m 






0.1 


0.340(4) 


0.3408(7) 


0.2 


0.485(5) 


0.4840(6) 


0.4 


0.693(5) 


0.6889(5) 



We can conclude from these data that our pre- 
dicted /ic moves with the quark mass in the ex- 
pected way; this behaviour increases our confi- 
dence on our results. In fig. 3 we report the chiral 



condensate at m=0.2 from our data and, super- 
imposed, the MDP results as reported in figure 2 
of Q. From fig. 3 it is evident a marked differ- 
ence between MDP results and those by our sim- 
ulations only for the largest lattice (8^ x 4) and 
again limited around Hc as in the SU (3) case. In 
particular the 8'^ x 4 results differ from ours at 
fj, = 0.6, the critical point derived in The au- 
thors of Q have tested the independence of their 
results on the initial configuration only for the 
4^ X 4 lattice. In this case MDP results agree 
at a good level with ours (see ipip of fig. 3). In 
our opinion the observed discrepancy has to be 
ascribed to convergence problems of the MDP al- 
gorithm, although they arise at volumes larger 
then in the SU{3) case. Once again there are 
serious doubts on the accuracy that the MDP al- 
gorithm can achieve near the critical region. 
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